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A spin- 1/2 chain model that includes three spin interactions can effectively describe the dynamics 
of two species of bosons trapped in an optical lattice with a triangular-ladder configuration. A 
perturbative theoretical approach and numerical study of its ground state is performed that reveals 
a rich variety of phases and criticalities. We identify phases with periodicity one, two or three, as 
well as critical points that belong in the same universality class as the Ising or the three-state Potts 
model. We establish a range of parameters, corresponding to a large degeneracy present between 
phases with period 2 and 3, that nests a gap less incommensurate chiral phase. 

PACS numbers: 05.30. Jp,75.30.Fv 



I. INTRODUCTION 

The proposal 1^ and subsequent realization of 
optical lattices for the manipulation of ultra-cold atoms 
has attracted a great deal of research towards the im- 
plementation of quantum computation ^jy,[3 and the 
simulation of condensed matter systems |3, S • 
The main advantages of optical lattices are their long 
decoherence times and high degrees of controllability. 
Hence, one is able to probe phenomena that manifest 
themselves at higher orders in perturbation theory such 
as many body interactions H^. This gives the possibil- 
ity to engineer and control exotic interactions in many 
body systems and subsequently realize novel phases of 
matter. Examples include cases where frustration effects 
are present or competing phenomena coexist giving rise 
to large degeneracy structures 0, . Towards this di- 
rection we consider a semi-one dimensional lattice com- 
prising a triangular ladder. The dynamics of two boson 
species mounting the optical lattice in the limit of strong 
collisional interactions can be effectively described by a 
chain of spin-1/2 interacting particles. Here we consider 
systems with two and three spin interactions given by 
= and ZZZ = Ei«+i<+2 respec- 

tively. We shall see explicitly that it is possible to make 
ZZZ dominant by appropriate tuning of the collisional 
and tunnelling couplings. As simple as these interactions 
may seem, when combined, they give rise to a rich variety 
of phase transitions, the study of which is the subject of 
this article. 

In particular, we consider a Hamiltonian that includes 
the two and three spin interactions with couplings Ai 
and A2 respectively, in the presence of a transverse field 
with unit amplitude. We observe that for A2 = the 



model is exactly solvable, as it reduces to the Ising in- 
teraction with critical points at Ai = ±1. For Ai = we 
obtain a Hamiltonian that is not analytically diagonaliz- 
able. Nevertheless, a numerical study reveals a trier it ical 
point at A2 = 1 that is in the same universality class as 
the three-state Potts model 16j. When the two spin in- 
teraction is dominant the ground state has a spin-order 
of period 2 while the dominance of the higher order in- 
teraction brings about a spin-order of period 3. Though 
these phases are gapped, there are values of relative cou- 
plings where these orders compete, giving rise to a high 
degree of degeneracy. This characteristic allows for the 
presence of a gapless incommensurate chiral phase that 
extends to a wide range of parameters. In Figure the 
fidelity of the actual ground state of the Hamiltonian is 
plotted against the ground state of each individual term 
comprising it. While we see that these agree for a large 
range of couplings there is an area between spin-order 2 
and 3 where these states fail to accurately describe the 
nature of the system. We shall identify this as the gap- 
less incommensurate chiral phase. A similar region has 
been presented by Fendley, Sengupta and Sachdev JJJ in 
a 1-dimensional hard boson system. 

The article is organized as follows. In Section ^ we 
present the realization of the two and three spin interac- 
tions by ultra-cold atoms superposed by optical lattices. 
Section Hill is a study of the tricritical point present at 
Ai = 0, A2 = 1. In Section HVl we present the properties 
of the full Hamiltonian. For that we employ perturbation 
theory, a bosonization of the model and a study of an ef- 
fective field theory that eventually reveal the incommen- 
surate chiral phase. Finally, in Section we conclude 
and present future directions. 
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FIG. 1: A plot of the fidelity, of the ground state \ip) 

with the predictions 1 0) obtained from each individual term of 
the Hamiltonian. The state 10) = |t|t ... |t|) corresponds 
to ZZ, 10) = llll ... til) corresponds to ZZZ and 10) = 
(ID — 11))'^'^ corresponds to the transverse field. The state 
10) — liii ••• iii) is the common ground state of ZZ and 
ZZZ for this range of couplings. The spin chain is of length 
N = 18. 

II. OPTICAL LATTICES AND THREE SPIN 
INTERACTIONS 

Consider the physical setup where an ultra-cold atomic 
cloud of two different species, a and 6, is superposed with 
a three dimensional optical lattice. The atoms can tunnel 
through the potential barriers of the lattice from one site 
to the next with a coupling J. When two or more atoms 
are present in the same site, they collide with coupling 
U. For sufficiently large intensities of the laser radiation, 
where J <^ the system is in the Mott insulator regime 
with a regular number of atoms per site. In particular, 
we can arrange the density of the atomic cloud to be low 
enough so that only one atom can exist at each site of the 
lattice, i.e. (n^ + n^) 1. In this way each lattice site is 
a simple two state system, that can be viewed as a spin- 
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FIG. 2: The one dimensional chain constructed out of equi- 
lateral triangles. Three-spin interaction terms appear, e.g. 
between sites i, i + 1 and i + 2 as for example tunnelling be- 
tween i and i + 2 can happen through two different paths, 
directly and through site i + the latter resulting into an 
exchange interaction between i and i + 2 that is influenced by 
the state of the site i + 1. 
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FIG. 3: The effective couplings Ai and A2 are plotted against 
J7[/ and J^/U for Uaa = Ubb = 2.12U and Uab = U. The 
coupUng Ai appears almost constant and zero as the unequal 
collisional couplings can create a plateau area for a wide range 
of the tunnelling couplings, while A2 can be varied freely to 
positive or negative values. 



1/2 particle. Within this representation, the Hamiltonian 
can be written in terms of Pauli spin operators. 

Indeed, for the triangular ladder seen in Figure [2l we 
obtain the Hamiltonian of the form 0, ^3 

H = B^< + Ai ^«i+A2 5]cTf<+i<2, (2.1) 

where 5, Ai and A2 are all functions of the initial tun- 
nelling and collisional couplings, J and U. The values 
of Ai and A2 can be controlled independently by varying 
J and U. In particular, it is possible to make A2 large 
compared to Ai as can be seen in Figure[2l A careful con- 
sideration of the system would reveal that terms of the 
form afaf_^2 ^^^o appear, due to the triangular configu- 
ration. This can be remedied by employing superlattices 
that do not alter the other terms. 



III. THE ZZZ INTERACTION 

Ultimately, we would like to study the entire (Ai,A2) 
plane, identify the regions of criticality behavior and un- 
derstand the global properties of the ground state. Before 
turning to the general problem let us consider the two 
special limiting cases. Initially, when A2 = 0, the Hamil- 
tonian reduces to the Ising interaction between neighbor- 
ing spins in the presence of a transverse magnetic field. 
This well-studied model exhibits criticality behavior for 
Ai = ±1. 

Another interesting model |0 can be obtained for 
Ai = 0. For simplicity we can take A2 = 1 and con- 
sider the Hamiltonian to be a function of the transverse 
field with amplitude B . As we vary B we observe 

that there are two distinctive regions. For B ^ 1^ the 
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ground state has all the spins oriented towards the x di- 
rection, and for 5 <C 1 there is a degeneracy between the 
states iTTiTTiTTi)^ its translations (the Z3 symmetry) 
and liiiiiiiii). It is of interest to study the behavior 
of the system between these two limiting cases. 

Due to the symmetry of the Hamiltonian we can 
pinpoint where a possible critical point can lie. For 



that, one can define, ji^ 



n 



n=0 



3n^ 



-3n-l- 



It is easily verified [2 



and 111 = 
that these 



operators obey the Pauli algebra commutation relations. 
Moreover, under this transformation the Hamiltonian be- 
comes BH[B~^) which has exactly the same spectrum 
as H{B). Hence, if there exists a single critical point it 
has to be at 5 = 1. 

We can verify this numerically and identify the critical 
exponents, by observing the minimum in the energy gap 
between the ground and first excited state on a finite 
chain of spins with periodic boundary conditions. In the 
thermodynamical limit this minimum, if it corresponds 
to a critical point, will become zero. Near this region the 
energy gap. A, is expected to scale as follows _22i] . 



(n^/^{B-Bc)). (3.1) 



Here <l> is a universal scaling function, z is the dynamic 
critical exponent, and v is the correlation length expo- 
nent. Figure El^a) gives a plot of Tyi oo^es^ versus B, 
which shows that the energy minimum progresses to- 
wards 5 = 1 for increasing A^. We may employ Eqn. 
(FrT|) to estimate z and v. Figure|Hb) shows Ari ooo2623^ 
against N^^^-^^^^^^[Bx — B^)^ where the data for systems 
of different sizes come together into a single curve. From 
Eqn. (|3.H) we deduce the value of the critical exponents 
to be z 1 and u ^ 0.76. 
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FIG. 4: (a) Quantum critical point signified by the minima 
of the energy gap A. Due to low these provide only an 
estimate for the position of the critical point. The progres- 
sion towards Ex = I 8iS N increases supports the theoretical 
predictions, (b) To determine the critical exponents we have 
chosen the values of z and u that best superpose the curves 
of Figure EJ a) into one uniform function 

We can also determine numerically the central charge 
of the critical theory and see if it corresponds to the same 
universality class. The central charge can be obtained by 
studying the scaling behavior of the entropy of entan- 
glement 22J. The latter is given by the von Neumann 



entropy, Sl, of the reduced density matrix, p^- It can be 
calculated from the density matrix of the original system, 
being in a pure state, where all but L contiguous spins 
are traced out. This indicates quantitatively the degree 
of entanglement of the L spins with the rest of the chain. 
Indeed, we have 



PL = tTN-L\i^)W, 

Sl = tr(pLlogpL), 



(3.2) 



where {ip) is the ground state of the system. For a critical 
configuration and for large L we expect Sl ^ log ^5 
where c is the central charge of the corresponding con- 
formal field theory and c is its complex conjugate. The 
central charge uniquely corresponds to the critical expo- 
nents of the energy and of the correlation length, z and 
i/, respectively. We know that for non-critical chains Sl 
should be saturated for large enough values of L. This 
behavior is observed from the simulations when A2 7^ 1. 
On the other hand, when we are at A2 = 5 = 1 we ob- 
tain Figure [51 that shows the expected logarithmic pro- 
gression. By a logarithmic fitting we can deduce that 
c ^ 4/5. Hence, our model is in the same universality 
class as the three-state Potts model and it corresponds 
to the critical exponents z = 1 and = 3/4 in agreement 
with the earlier findings. 
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FIG. 5: Entropy of entanglement, S*!,, as a function of L for 
a total of 19 spins. The plot shows a logarithmic behavior 
that indicates criticality for A2 = 1. 

IV. THE FULL HAMILTONIAN 

We can now turn to the full Hamiltonian for arbitrary 
values of the coupling parameters. For convenience we 
set B = 1 and vary only the interaction couplings Ai 
and A2. Without loss of generality we can restrict on the 
A2 > half-plane as the case A2 < is automatically 
obtained by exchanging ^< — >l. 

For large values of Ai or A2 it is possible to neglect 
the transverse field X = ^ • af and treat the system 
classically. We find two asymptotes along which the 
ground state should undergo a first order phase transi- 
tion. By introducing X as a small perturbation one can 
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predict how these curves behave when we approach the 
origin Ai = A2 = 0. We find that along the asymptote 
A1/A2 = 3/2 the classical ground state is highly degen- 
erate: it is in fact infinitely degenerate in the thermo- 
dynamic limit. In the quantum regime, this gives rise 
to an intermediate region between the two phases with 
spin-order 2 and 3, the size and nature of which will be 
studied in this section. 



A. Perturbation Theory 

To understand the behavior of Hamiltonian (|2.1|) we 
shall employ perturbation theory. Surprisingly, the sec- 
ond order perturbation will give a very good approxima- 
tion to the numerical findings shown in Figure Q 



1. Classical Regime 

As a first step we take Ai, A2 ^ 1, thus neglecting the 
X term the Hamiltonian reduces to 



H ^ Ai ^crfcrf+i + A2 



(4.1) 



This Hamiltonian has a classical behavior, so it can be 
solved exactly by examining eigenstates of crf^. There 
are three regions on the (Ai, A2) plane of distinct ground 
state behavior. Consider first Ai < 0, A2 > 0. The energy 
of the system is minimized uniquely by the state 



|A> = Uiiiii..4). 



(4.2) 



From the expectation values, {A\'^-a^a^_^i\A)/N = 1 
and {A\'^-a?af_^ia?_^2\^)/^ — can calculate 

the energy per site to be Ea = Ai — A2. Similarly we 
take the cases of Ai > 0, A2 ^ Ai and Ai > 0, A2 <C Ai, 
resulting in the following ground states 



|5) = lTnTn...i), 

with Eb = — 1^1 — ^2 and 

\c) = \nnn---i), 



(4.3) 



(4.4) 



with Ec = — Ai. Observe that these states are 3 and 
2 periodic respectively giving rise to Z3 and Z2 symme- 
tries. To evaluate the position of the phase transitions 
we simply equate energies in neighboring regions. In this 
way we obtain the asymptotes Ai = and Ai = |A2 that 
separate the states of the system \B) and \C) as seen 
in Figure El 

Another classical region is obtained at Ai,A2 <C 1, 
where the transverse field, X, dominates. In this region 
the ground state is given by 



\D) = 



(4.5) 



with |— ) = (It) — \ i))/V^ and an energy per site of E^ 



2. Quantum Regime 

We now consider the transverse field as a perturbation 
to the theory, so we rewrite our Hamiltonian as 



H = X' 



+ 1^2 + 2 



), (4.6) 



where A = I/A2 and fi = A1/A2. Cyclic permutations of 
states \B) and \C) are also ground states of the classical 
theory due to the translational invariance of the Hamil- 
tonian. Nevertheless, we need not employ degenerate 
perturbation theory at this point since the permutations 
behave identically under the transverse field, X. For this 
Hamiltonian the energies per site are 



Ea = Ai 



Ao 



1 1 

A^6-4/i' 



Eb 



Ai + ^)-A2 



3 J 6 + 4/i 



(4.7) 



. (4., 



By inspection we see that the two energies are equal for 
/i = 0, which coincides exactly with the classical solu- 
tion. Hence, at 0(A^), the boundary between \A) and 
\B) remains unchanged, as seen in Figure El 

Let us turn to the perturbation treatment of the 
boundary between the \D) and the \A) states. In fact, to 
the second order in the couplings, Ai and A2, we find 



E 



D 



A 



(4.9) 



-1. 



Equating Ejj = Ea we obtain the corresponding bound- 
ary to be a correction of the classical one given by 
Ai = A2 — 1. The computational solution of Eb = Ea 
is given in Figure El It provides a line of second order 
phase transition joining the Potts critical point and the 
Ising critical point that fits firmly with the fidelity plot 
in Figure n 



3. Competing Ground States 

We saw above that for /i = | the period 2 and period 3 
phases have the same ground state energy. Interestingly, 
on this asymptote the two spin-orders can mix, giving 
rise in the thermodynamic limit to an infinitely degener- 
ate ground state. A perturbation theory around this area 
is non-trivial as care has to be taken in the way the de- 
generacy is lifted. Using a bosonization procedure similar 
to that presented in 14], we will calculate the spectrum 
of the theory up to second order in the transverse field. 
We will find a finite intermediate phase separating the 
period 2 and period 3 gapped phases. 

Let us present in detail the bosonization procedure nec- 
essary to develop the perturbation theory. Consider two 
new types of bosons - the |3), created by the operator. 
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and 1 2), created by . In the spin formalism, |3) cor- 
responds to iTTi), while 1 2) is equivalent to |||). Thus, 



\B) = I33...3) and \C) = |22...2). 



(4.10) 



All states that are composed of 2s and 3s are degenerate 
along the line /i = §• Therefore, the region of the phase 
diagram that we want to study is given by 



0<|/i-3/2|<|Ai|,|A2| 



(4.11) 



Since there are no first order contributions in the per- 
turbation theory, the region we want to study can be 
parameterized by 



(4.12) 



for a dimensionless parameter a. 

The approach is to create a theory with t and p bosons 
that is equivalent to the present one. Let us first work 
with t bosons. The vacuum state has the form |22 . . . 2), 
with no 3s. To construct the Hamiltonian we need the 
self energy of a t boson, the interaction energy and the 
hopping amplitude. For that, we calculate the energies 
per site needed to add 22, 33 and 23 bonds, given by E22, 
Ess and E'23, respectively. In particular. 



I 22... 2; 



m—l 



1 22... 22; 



(4.13) 



gives rise to the energy gap £"22, where m is an integer. 
A similar transformation gives E^^. Furthermore, 



. 2323 ) 



(4.14) 



induces the energy gap 2£^23- The energy of each of 
these states can be calculated perturbatively. Up to order 
0(A^), we obtain 



22 



-2/i-h 



2(1 -2/i)' 



E: 



33 



^23 



X2 

-(3 + /i)-y, 

-2(^+^)-y 



(4.15) 



Consider now the energy, Et, of creating a single t boson 
in a background of 2s, and the interaction energy, Et^int: 
generated when two t bosons are brought to adjacent 
sites. The former can be evaluated by considering the 
energy gap of the transformation 



I 22... 222 i 



5m 



1 23. ..23), 



(4.16) 
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which creates 2m bosons. The states are chosen so that 
the lengths of the underlying spin chains remain the 



same. For the interaction energy, the relevant energy 
gap arises from 



|23_^^) ^ |2_^3_^). 

4m 2m 2m 

Hence these energies per boson are given by 



(4.17) 



Et 



2E2S — -E22 



+ 2a 



E, 



't,int 



2E' 



23 



Eo 



E: 



33 



a2 
4 ■ 



(4.18) 



Finally, we calculate the hopping amplitude, which is 
given by the activation energy needed to perform the 
exchange 

|23) ^ |32). 

This has a value of — We are now in a position to 
rewrite the Hamiltonian, which takes the following form 

A2 



-2/i 



2(1 - 2/i) ^ 



(4.19) 



24 



+ 2cr + 



1 t t 



The label / is over the original indices, with a t boson 
being centered in the middle of the 3 deformation, and 
the integer is the total number of sites in the original 
spin formalism. We consider cr ^ 1, which corresponds 
to the spin-order 2 phase. The ground state is the vac- 
uum state |22 . . . 2), which is in agreement with original 
analysis. By employing the Fourier transform of ti we 
obtain that the lowest excited energy above the vacuum 
state is A^ (2(j — Thus at cr = ^ a phase transition 
occurs from period 2 to the intermediate region. 

Following the same procedure we can find the bound- 
ary of the period 3 gapped phase. Now the vacuum state 
is |33 . . . 3) and the relevant energies are with respect to 
the p boson. We have 



En 



2E23 — -Ess 



A2 



E, 



'p,int 



2E2S — E22 — Ess 



4 2 

a2 



(4.20) 



while the hopping amphtude is the same. In terms of 
these variables the Hamiltonian takes the form 



H 



N A' 



->^'Y.\\{pI^P^+p\p^+^) (4.21) 



-^Pi+2PiPl+2Pl 
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in view of Eqn. (|3.1|) . we confirm that as the transverse 
field becomes insignificant, the phase transitions return 
to first order behavior with an actual energy level cross- 
ing. Along the upper asymptote, Li, we observe similar 
behavior. The critical exponent is large in configurations 
far away from the origin and decreases continuously down 
to z = 1/2 at the Ising critical point Pi. 

These results present a rich variety of phase transi- 
tion behavior. From Figures Q and El we see that the 
second order perturbation theory gives a good picture of 
the ground state of the full Hamiltonian in most of the 
parametric region. The only part of the (Ai,A2) plane 
that is not yet understood is the strip bounded by the 
lines, Li and L2. Probing the physics of this region is 
the subject of the next subsection. 



FIG. 6: An outline of the perturbative results is given by 
the solid lines. Dotted lines are inserted by hand for comple- 
tion. There are three disjoint regions of different periodicity, 
separated by lines of phase transitions, of both first (L3) and 
second order (Li, L2 and L4). 



An identical analysis is now necessary, by simply con- 
sidering (J <C —1. The p bosons cost large positive 
energy, making the ground state the vacuum, as one 
would expect. The first excited state has an energy of 
— + 1), so the boundary is at a = — 

We therefore have a prediction of the width of the in- 
termediate phase that is illustrated in Figure This 
result is in perfect agreement with the fidelity plot in 
Figure 



C. Incommensurate Phase 

It becomes apparent from the previous that perturba- 
tion theory cannot give reliable information about the 
physics between Li and L2. To study that region we 
resort to the Landau- Ginzburg approximation As 
we shall see, the resulting theory corresponds to the chi- 
ral clock model [2^ , which predicts non-zero chirality for 
our model. This also indicates the presence of a gap- 
less incommensurate phase between Li and L2. For a 
finite chain we find strong evidence for the presence of 
this phase using exact numerical diagonalization. 

To gain a further insight into the nature of the ground 
states of our model we introduce the order parameter Q 



li),UI, 



(4.22) 



B. Criticality Behavior 

At this point we can discuss the criticality behavior 
of our system. Apart from the Ising critical points. Pi 
and P2, and the three-state Potts model critical point, 
P3, we can identify the critical behavior of the lines, Li, 
L2, L3 and L4, seen in Figure El where our results of the 
perturbation treatment are summarized. 

Specifically, we obtain first order critical behavior for 
Ai = and A2 > 1 and second order behaviors for the 
rest of the curves. The perturbation theory shows that 
the line, L3, acquires no second order contribution indi- 
cating its first order nature. This is also supported by a 
numerical study of the energy gap for a chain of 18 spins. 
Indeed, we find that its minima lie along lines similar to 
the ones in Figure El Furthermore, there is an actual en- 
ergy crossing along the first order line while all the other 
critical points remain gapped. 

Studying the critical exponents along the line L4 we 
find that in moving from P3 to P2 along the line of crit- 
icality, the exponent z decreases to zero and becomes 
discontinuous at the point P2. Moving in the other di- 
rection along L2 yields ever increasing values of z. Hence, 



where the operator is the projector onto the down 

state of the j*^ site. The periodicity of a state is re- 
vealed by the expectation value of it is maximal for 
p-periodic states, and minimal (or zero) otherwise. In 
particular, the wave order parameter, ^3, is a complex 
operator that detects the period 3 states. We wish to 
use ^3 as an order parameter to probe the behavior of 
our model between the boundaries Li and L2. With this 
in mind we construct a continuum quantum field theory 
with an action that has the same symmetries as (^3) 
does in the Z3 region. 

It is easy to verify that (^3) is invariant under trans- 
lations ^3 e^^*^/^^3, for integer I. A more intriguing 
symmetry is given when simultaneous spatial refiections 
and complex conjugation are performed, x — x, ^3 
^3 Finally, (^3) is invariant under time reversal, 

r — r. According to these symmetries, we can write 
down the corresponding effective field theory, given by, 



m^^S^^3 + c.c. + |S^^3 



Ss = J dxdr 



r 1^3 



c.c. 



(4.23) 
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FIG. 7: Surface plot of (^18/7) • There is a narrow band of 
excitation within the intermediate region predicted by per- 
turbation theory, giving strong indication for the presence of 
gapless incommensurate phase. 

The parameters a, v and r are non-universal functions of 
the original Ai and A2. It has been shown 27] that this 
model has a critical point at a = r = that belongs 
in the same universality class as the point P3. We may 
therefore identify them. 

The term a^!%^dx^^ breaks chiral symmetry when ^3 
is complex. The resulting theory has been studied by 
Ostlund and it has been shown that it belongs in the 
same universality class as the chiral clock model. Hence, 
it suggests the presence of a chiral phase for non-zero 
values of a. Furthermore this model predicts that this 
phase is gapless and incommensurate 26, 28]. By anal- 
ogy we conjecture that our model has such a phase away 
from P3 and we shall investigate it numerically. 

The gapless and incommensurate characteristics of this 
phase indicate that (^p) can build a non-zero value by 
applying the appropriate perturbation. Hence, we add to 
the Hamiltonian a small periodic potential of the form 

V = lQ-^Y,e^^'i'^\i).{i\ (4.24) 

3 

and we plot (^p) for several p as a function of Ai and A2. 
Indeed, for p < 3, we find that (^p) 7^ between Li and 
L2 and deduce that a 7^ there (see Figure [7|. A naive 
analysis of our field theory in which we treat as a 
fluctuation of the uniform solution indicates that a may 
be proportional to 27r/p — 27r/3. This suggests a region 
of a 7^ 0, with opposite sign corresponding to p > 3, in 
which a second chiral phase exists on the other side of 
P3. We discover numerical evidence for this in Figure [HI 
when Ai < and A2 < 1, but its effect is of significantly 
lower order than in Figure [7| 

Along with the gapless incommensurate properties of 
this phase we would like to investigate its chiral nature. 
Note that, on the one hand, the chirality operator, given 
by Xi = ' <^i+i X is an imaginary hermitian op- 

erator. On the other hand the Hamiltonian (|2.1|) is real. 



FIG. 8: This plot shows (^is/s - ^18/7) over Ai < 0,A2 < 
1. There is a clear region, bounded by the predictions of 
perturbation theory, where responds best to excitations 
of spin-order p = 18/5 > 3, indicating, in tandem with Figure 
a theory with p < 3 above P3 and p > 3 below. 



Thus, if the Hamiltonian possesses a ground state with 
non-zero expectation value (x) , that ground state should 
be mult idegener ate. If this is the case in the thermo- 
dynamic limit we would like to see if it is also approx- 
imated in the finite case. By numerical diagonalization 
of the finite-size Hamiltonian, we can obtain the lowest 
energy gap, AE^ for values of the couplings Ai and A2 
between the Li and L2 lines. In Figure [HI we can clearly 
observe an exponential damping of AE as a function of 
the length of the spin chain. This allows us to assume 
that in the thermodynamic limit the spectrum will even- 
tually become degenerate. 

0.06^ 




I ^ ^ ^ ^ ^ 1 

8 10 12 14 16 18 20 
N 

FIG. 9: The energy gap, AE, per site between the ground 
state and first excited state tends exponentially to zero as the 
total length of the chain increases. In the thermodynamical 
limit we expect to have zero energy gap. The lines correspond 
to different points along the A1/A2 = 3/2 asymptote, ordered 
from top to bottom the further away the points are from the 
origin. 
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V. CONCLUSION 

The present spin model approximates the behavior of 
the Mott insulator of two species of atoms when a tri- 
angular ladder geometry is realized. The spin interpre- 
tation, which makes the study of the model simpler and 
more intuitive, is valid when there is only one atom per 
lattice site, i.e. at the regime of strong collisional cou- 
plings. With this condition satisfied the criticality be- 
havior of the spin model straightforwardly represents the 
behavior of the Mott insulator. We can argue that non- 
zero chirality in the spin model corresponds to a ground 
state with persistent currents . This suggests a coun- 
terfiow of the two different atomic species that could be 
measured by using atomic spatial correlations j3G] . 

To summarize, in this article, we have generalized the 
one-dimensional Ising model, ZZ, to include an addi- 



tional triple interaction, ZZZ, in the presence of a 
transverse magnetic field. A rich criticality behavior has 
been revealed when the coupling of ZZ and ZZZ are var- 
ied. The phases of the ground state are identified accord- 
ing to their periodic structure. In particular, a complex 
order parameter related to the states with period three 
has revealed a chiral, gapped, incommensurate phase. In 
[31 ^ thorough study of this phase in a one-dimensional 
hard boson model exhibiting similar behavior was per- 
formed. An in depth study of the incommensurate phase 
will be the subject of a future publication. 
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